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Abstract 

We compare two non-relativistic (NR) reduction schemes (heavy-fermion and Foldy-Wouthuysen) 
that are used to derive low-energy effective-field-theory Lagrangians. We give the explicit trans- 
formation between the two types of fields to ©(l/m^), derived from a quite general, relativistic 
Lagrangian. Beyond leading order the NR reductions always involve the smaller components of the 
Dirac spinors that are to be integrated out to formulate the NR theory. Even so, the transformation 
between the NR Lagrangians can be carried out explicitly to 0(l/m^) using a field renormalization, 
as long as the lower components of the Lagrangian are known. The fixed coefficient corrections 
to some low-energy constants at 0(l/m^) will depend on the particular scheme chosen, but will 
match after the field renormalization. 
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I. INTRODUCTION 



When considering low-energy phenomena, it often simphfies the calculations to start from 
a non-relativistic Lagrangian. Several useful and well-known examples can be found. These 
include the low-energy few-body nuclear systems for which chiral perturbation theory (xPT) 
has been developed, the electromagnetic properties of composite particles (moments of the 
proton and neutron), atomic physics [non-relativistic quantum electrodynamics (NRQED)] 
and heavy-quark effective field theory (HQEFT). Since the original theory, here assumed to 
be an effective field theory (EFT) in the Weinberg sense is generally given by a relativistic 
Lagrangian, one faces the problem of performing a non-relativistic (NR) reduction. 

Several solutions have been suggested over the years. The simplest, which we will label 
the direct Pauli reduction scheme or the 1/m expansion, expands the Dirac matrix element 
in terms of the two-dimensional Pauli-spinors that appear in the usual representation of the 
four-spinors. 

HQEFT was originally formulated in the non-relativistic limit, utilizing the constraints 
given by the symmetries of quantum chromodynamics (QCD) in the static limit 2|. Later 
HQEFT was elegantly derived from the relativistic QCD Lagrangian using a path-integral 
method j^; this we will call the heavy-fermion (HF) approach (for a recent review of HQEFT, 
see jj]). These ideas were carried over to the descriptions of heavy baryon chiral perturbation 
theory (HB^PT) 

An alternative NR reduction method would be the classic Foldy-Wouthuysen (FW) trans- 



formation 



B3, 



which was designed to provide the NR limit of quantum electrodynamics 



(QED). It has also been used to derive HQEFT 10|, but this work was initially criticized by 
'or not proving HQEFT to be an EFT (see, however, later work 111, [l2])- According 



Ref. [3| 

to Ref. [3j, HF and FW give identical results for HQEFT, although as far as we are able 
to ascertain, this was shown only to (9 (1/m), where m is the (heavy) fermion mass. Later, 
these two methods were indeed shown to be equivalent, after a wave function renormaliza- 
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16[}. However, these references 



tion and a unitary transformation Q, [l^ (see also Refs. 
considered only renormalizable theories (QED and QCD) and 1J| only the explicit example 
of a constant electric field (in QED). Neither considered the low-energy constants (LECs), 
which parametrize the (usually unknown) short-distance physics in an EFT like, e.g., xPT. 
It has already been shown that, when the relativistic theory is fully known, the direct 



2 



Pauli reduction scheme (1/m expansion) differs from FW and only the latter gives S'-matrices 
that are in agreement with the relativistic ones [l^]. This provided part of the inspiration 
for the present investigation. Another motivation was the claims in the literature that FW 
and HF are equivalent for HQEFT |3|,Q as well as for HBxPT Q. 

In both HF and FW the contributions involving the lower Dirac spinor components are 
eliminated, and this can be formalized by integrating out these "small" components from 
the path integral. For momenta and field energies smaller than m, the surviving terms can 
be expanded in a perturbative series. The difference between the two methods lies in the 
way the separation of lower and upper components is made, resulting in different NR fields. 
Note that this 'integrating-out' procedure is different from what is done in, e.g., electroweak 
theory, where the heavy vector bosons are integrated out completely, at lowest order giving 
the Fermi weak interaction model. Here we rather want to keep the heavy (fermion) degrees 
of freedom (when they are protected by a symmetry, like baryon number conservation) , but 
eliminate the explicit antifermion contributions (lower components), since they can not be 
excited in the non-relativistic limit. 

In this paper we will study the HF and FW representations in detail, starting from a quite 
general, relativistic, Lagrangian, which is assumed to be known. The two methods appear to 
disagree at higher orders in 1/m, i.e., the so-called "fixed-coefficient" terms (1/m corrections) 
differ between the Lagrangians. We will give explicit expressions for the transformation 
necessary to get from HF to FW to (9(l/m^), showing that the corresponding fields each 
contain pieces of the small/lower components of the other, already at 0{l/m). Even so, the 
functional form of the corresponding Lagrangians will differ first at (9(l/m^), which explains 
the equivalence found between the 0{l/m) results of Refs. jsl and 10 1. It is possible to 
show the equivalence of the two approaches [l^] (at least for a constant electric field in 
NRQED), by using field renormalization and unitary transformations. We want to treat a 
more general Lagrangian than what was done in |l^ and will show that at higher orders there 
are differences in the fixed-coefficient (1/m) corrections to the LECs in, e.g., xPT. (In this 
paper we will only be concerned with LECs introduced already in the relativistic Lagrangian, 



as done in, e.g.. 
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20|.) Thus, one has to be careful to use correctly renormalized fields 



when comparing LECs at this order. At higher order also unitary transformations need to 
be applied. 

This paper is organized as follows: In Sec. [Ill we give a brief review of the essential 
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features of the HF and FW schemes, followed by a detailed comparison and derivation of 
the transformation between them. The paper is concluded by a discussion in Sec. IIIII Some 
mathematical details are collected in an Appendix. 



II. NON-RELATIVISTIC REDUCTIONS 



We will assume as the starting point the relativistic Lagrangian 



C = ij{i'p-m + 7oG) ip, 



(1) 



where the covariant derivative "D^ = + and the general coupling matrix G (the extra 
factor of 7o is included for future convenience) is 



where A-D are 2x2 matrices. The connection F'^ and the coupling matrix G contain all 
single-fermion couplings (to, e.g., pions and electroweak fields) allowed by relevant symme- 
tries, e.g., chiral symmetry, Lorentz invariance, parity, flavor and color SU(3) etc. Through- 
out the paper we will work in the usual Dirac representation and G should be interpreted 
in that representation. We further assume that G is hermitian, i.e., that A and D are her- 
mitian and B'^ = G. No assumptions will be made regarding the commutative properties 
of the elements of G and V, so the results of this paper apply equally well to Abelian and 
non-Abelian fleld theories. Note that A-D and will in general contain higher order 
interactions terms which should be expanded in the final NR formulas. This expansion is 
separate from the 1/m expansion created by the NR reduction, but in some cases the two 
expansions can be combined into a single counting scheme. This happens for, e.g., HB^FT, 
where m ~ A^ and A^ ~ 1 GeV is the chiral scale. The generating functional for the 
Lagrangian ([1]) is 



where X represents the fields implicit in F^ and G. These fields will from now on be 
suppressed in our notation. 




(2) 




(3) 
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A. A brief overview of the heavy-fermion method 



If the fermion is heavy, we can for low four-momenta treat it as essentially static and 
expand around its large mass. Following Ref . , this is accomplished by writing the fermion 
momentum p^^ in terms of the four-velocity f ^ and a residual fermion momentum l'^ as 

p^' = mv^ + (4) 

where <^m and = \. Consequently, for an on-shell fermion, = implies that 

2mv -1 + ^ = 0. (5) 

The fermion field is next split into large (H) and small (h) components 

iIj = e-'"'''''{H + h), (6) 

where i/H = H and t/h = —h, which shifts all the linear mass dependence to the small 
component. Without loss of generality we will assume for the remainder of this article that 
v'^ = (1, 0, 0, 0), and thus t/ = 70, which simplifies the algebra and the comparison with FW. 



(It is possible to keep general and also do the FW transformation in terms of 10|, ll3| , 
but we choose a simpler approach here.) The resulting Lagrangian is 

^ ^ y-ia ■ V + C 2m + iVo + D J y h J 

where we have chosen to work with instead of H to ease the comparison with the 
FW expressions further on. The cross talk between upper and lower components can be 
eliminated by defining a new small component field h' Sj such that 

h = h' -{2m + iVo + D)-\-ia ■ V + C)H, 

= h'^ - {-ia-V + B) {2m + Wo + D)-\ (8) 

This turns the Lagrangian into a block-diagonal form: 

C = H^ (iVo + A- (~ia ■ V + B) \- -{-ia ■V + C)\h + h'\2m + iVo + D)h', 

\ 2m + iVq + D J 

(9) 

and 'completes the square' in the path integral exponential. Up to this point no approxima- 
tions have been made, so this equation can be regarded as a Dirac equation in two-component 



form. After 'integrating out' the small components, which gives a constant determinant [3| 
(see the Appendix) , only the H fields remain and we have obtained a rudimentary form of the 
HF Lagrangian. Note that the H fields have to be renormalized to preserve the norm 
The antifermion is implicitly included as the '2;-graph'-like term of Eq. We next expand 
this term, assuming that iD^^ G <^ 2m. Thus, this part of the HF Lagrangian is expanded 
in powers of {iT>Q + D)/2m, which contains only time derivatives of the fermion field (apart 
from possible derivatives in G), sandwiched between two space derivatives. The first few 
orders of the HF Lagrangian are then 

= H\iVo + A)H, (10) 
= -l-H\-ia-V + B){-ia-V + G)H, (11) 

= A_H^^-ta-V + B){Wo + D){-ia-V + C)H, (12) 

where the indices indicate inverse powers of m. The Lagrangian in general also contains 
higher order interaction terms in A-D and V^. This expansion will not be shown explicitly. 
For the purposes of this paper we count X^q ~ ~ 111; although counting Vq ~ F/m is also 
possible, depending on the process being studied. Short-range interactions are parametrized 
as local counter terms, which contain LECs. The LECs are determined either by matching 
conditions (when the underlying, more complete theory is well known, e.g., HQEFT and non- 
relativistic QED (NRQED) [22]) or by fitting EFT calculations to data (when the underlying 
theory is not well known, e.g., HB^PT). In some cases, models such as resonance saturation 
are used to evaluate or estimate the LECs. 

The propagator for the H field is given by Eq. (Q 

' + (13) 



'-' 2m+ido ^ 2m 

where □ = 9^ - The inverse of this propagator is non-linear in ido beyond a certain 
order, and this propagator usually appears in expanded form . In order to compare 

this HF propagator to the FW propagator that we will give below, we choose a situation that 
appears often in our own line of research, i.e., pion-nucleon dynamics in HB^PT. Consider 
an on-shell nucleon with residual four-momentum /. The HF nucleon propagator after it has 
emitted one pion of four-momentum q (real or virtual, but with qo, |q| <^ m) can be written 
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as 



(1+^) _ ^(1 + ^) 



7 n ^ (^-q)" nn - + 9' ' 

- go + -Qo ^ + 2^ 

where the relativistic on-shell condition ([5]) has been apphed. The propagator f|T^ can be 
further expanded as needed, whether in 1/m or in ^JrrL^^Jm as in pion production off two 



nucleons M(^h,s application of the propagator rs smrUa.^ to the orre advocated by Harrhart 
and Wirzba [24]. After a simple re-identification of I and q as the electron and photon four- 
momenta, Eq. (IT^ can also represent the NRQED electron propagator after emission of a 
soft photon. 

The HF propagator [Eq. ( fT4l) ] is identical to an 1/m expansion of the relativistic propa- 
gator, projected on the upper left 2x2 submatrix: 

iimil + / — ^ + va) . Im + ~ Q'o 



[mv + I — qy — m'^ —2{m + lo)qo + 21 ■ q + 



1 _|_ Iq—Qo 
^ 2m 



n ^'1 4- 1 

-Qo - — + 2^ 



(15) 



B. Foldy-Wouthuysen transformation in a Lagrangian formulation 

In order to compare to the HF results we choose to derive the FW transformation from 
the Lagrangian |lO, [l^ instead of the Hamiltonian as is usually done. In addition this 
approach makes it possible to treat the time and space components in a more symmetrical 
manner. The time derivative is implicit in all expressions containing Vq. Our starting point 
is the Lagrangian ([T]) rewritten as 

£ = (^q^^ _ .v_^^rn + G^ (16) 

where the field tp"^ is chosen instead of since it makes the algebra somewhat more trans- 
parent for the present purposes. This expression can be further separated into odd and even 
terms, i.e., terms that mix and do not mix the upper and lower components (in the Dirac 
representation) 

C = i;\-S-0-^om)^, (17) 

where the sign convention for the even (S) and odd (O) operators is in accordance with 
Bjorken-Drell j^. The mass term can be removed from the upper components by shifting 
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the zero-point of the energy, i.e., we define if) = e giving 

C = ^^[-8 -O - (70 -l)m]ij. (18) 
The exphcit expressions for the even and odd operators are 



£ - - . (19) 




O ^ - _ _ . (20) 



The FW transformation redefines the fermion field through the consecutive apphcation 
of unitary matrices: 

= e''"e'''e'^i,, (21) 

where 5" is related to O' in the same way as 5* to O [Eq. fl22l) ] and similarly for S" and O" . 
These higher order odd operators are introduced below. The transformations carried out to 
0{l/m?) lead to the Lagrangian 

(23) 
(24) 



C" 


= i,'"^ {-8' - 


(7o-l)m)*l"'. 


8' 






O' 


= ^[0,£]- 

2m 




O" 


2m 





(25) 
(26) 

where it is assumed that the time derivative in 8, i.e., in Vq, acts on everything to its right. 
All odd terms have been transformed away to this order and only the even operators remain. 
The explicit expressions for the higher order odd operators are 

if B'\ 
2m J 

^ ( B'iWo + D)- {^Vo + A)B' 

4m2 1 + A)- {Wo + D)C' 
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where we have defined 



B' = {-la ■ V + B){iVq + D)- {Wo + A){-ia -V + B), (29) 
C' = {-ta ■ V + C){iVo + A)- [iVo + D){-ia -V + C). (30) 

The upper (A^) and lower {n) FW components of if)'" [Eq. f|2ip ] are related to the HF fields 
H and h by the transformation 

1 — TT^BC + TT^B'C — Tr~-B + T^B' + ,^ J BCB — -^^B" 

8m^ 8 m'' 2 m 4m^ 16m'' 8m'' 

1 Q _j_ 1 ^/ _ 3 (J Q(J -|- 1 (J" ]^ —}—(JQ _ —}—(J'B 

2m 4m^ 16 m'' Sm*' 8m^ Sm'' 







I „ J 





where we use the abbreviated notation 

B = -ia-V + B, (32) 

C = -ia-V + C, (33) 

B" = B'iiVo + D)- {Wo + A)B', (34) 

C" = C'iiVo + A)- {Wo + D)C'. (35) 

The Lagrangian after the three FW transformations is 

dfv^ ~ 2mn^n, (36) 

4w = ^^(^^0 + A)N + n^{iVo + D)n, (37) 
4w = -^N\-ia-V + B){-ia-V + C)N 

+ -^n\-ia -V + C){-ia -V + B)n, (38) 



4'w = -^N\-ia ■V + B){iVo + D){-ia ■V + C)N 

- -^N^{{-iB ■ V + B){-%d -V + C), {iVo + A)}N 
+ -:^n^ {-i^ ■ V + C){Wo + A){-ia ■ V + B)n 

- ^n^{{-i^-'D + C){-ia-V + B),{iVo + D)}n. (39) 

By construction terms mixing upper and lower components start to appear only at the 
next higher order, so the square is already 'completed' and one can integrate out the rt^n 
contributions to this order. As in the HF case, the integration over the small components will 



give a constant determinant (see the Appendix). When higher order terms, involving higher 
order LECs, are included in the relativistic Lagrangian, one has to take into account that 
A-D and contain higher order interaction terms, such that further expansion of these 
quantities becomes necessary. Just like HF, the FW transformation generates an expansion 
in 1/m, but the two series are not identical. 

For example, the free N propagator is (see Ref. 

(40) 



i.e., exactly what would naively be expected for a non-relativistic propagator. Unlike the 
inverse of the HF propagator f[T51) the inverse of fHUl) is linear in ido, although its dependence 
on the spatial components is more complicated. Applying this FW propagator to the same 
situation as in the HF case (but with the relativistic on-shell condition now being Iq = 
V^r^^TV -m= ^- ^...) we get 

i = (41) 

h-qo 2^ + + • • • -QO 2^ ^ 8^ ^ • • • 

which differs from the HF propagator [Eq. (1141) ]. Although the HF and FW propagators 
go on-shell for the same value of /'^ for any given order, in general they do differ off-shell 
beyond leading order. Also the on-shell residues are different, as discussed already in [l4|. 
It was pointed out in Ref. |l4], however, that the difference between these propagators can 
be compensated for by a field redefinition. 



C. Comparison of HF and FW to third order 

In contrast to HF, FW is an expansion in terms of both space and time derivatives of 
the fermion, i.e., the FW transformation is not restricted to the time derivatives only. The 
two schemes also differ in that the integrated-out pieces of FW have contributions at all 
orders, not just the first few as in HF — compare Eq. iQ to Eqs. fl36l) - fl39l) . Comparison with 
Eq. (IT^ shows that the difference in the Lagrangian between these two expansion choices 
is given (to this order) by the second line of Eq. fl5^ . i.e., differences start to appear at 
0(l/m2). Both result in expansions in 1/m, where m indicates the size of the cutoff of 
the theory. Note that the fields H and are not identical (and neither are h' and n), 
except at the lowest order, but are related through a rather involved transformation, see 

10 



below. In contrast to the transformation in [l^, which was carried out to 0{l/nT') for a 
constant electric field, we have kept the Lagrangian general (not necessarily renormalizable, 
e.g., HBxPT), but restricted the expansion to 0{l/rn?). 

In order to transform from the HF to FW formalisms we need to 'complete the square' 
for the H and h fields according to Eq. ([H]). This boils down to an additional transformation 

(42) 





H 


\h j 





2m+iT>o+D 

We then arrive at the field transformation 




n 




^/-W = I I =M« I , (43) 



where 




-^B{W, + D)C ihBCB - ^B" 

-^CBC + ^[C - {^V, + D)C]{^Vo + A) -^C'B 



(46) 



We have checked that this transformation indeed transforms the Lagrangian from the FW 
to the HF representation [to (9(l/m^)]. This is not a unitary transformation, because of the 
block-diagonalizing done in Eq. (1421) . It is equivalent to the wave function renormalization 
and unitary field redefinitions in Refs. |2l|. In contrast to Ref. [14], however, we give 



the general transformation to (9(l/m^), i.e., it applies equally well to non-Abelian and non- 
renormalizable theories. Obviously, the FW spinor contains parts of the (to be) integrated- 
out HF spinor h'. This means that transforming the upper/large components from one to 
the other mixes in lower components of the other representation. Thus, the transformation 
can be carried out explicitly only in cases where the lower component Lagrangian is known 
(up to the required order). 

The general relation between the HF and FW Lagrangians to a given order is 

H^UhhH + h'^L^^H'h' = N^UmnN + n^L^nU. (47) 

11 



Assuming a general form for the ith order matrix M^"^' = that relates the FW 

and HF fields [Eq. P3|) ] we can write [hereafter suppressing the index (z)] 




at ct 




6t c/t 

This puts constraints on the Lagrangian operators Uxx and L^x- 



(48) 



Uhh = a^UNNa + c'LnnC, (49) 

-^^h'/i' = b^UNNb + d^Lnnd, (50) 

= a'UNNb + Lnnd, (51) 

= b^UNNa + d^LnnC. (52) 



By taking advantage of the block-diagonal conditions [(ETj) and ( ]52l) ]. the lower component 
L„„ can be eliminated from the expression for Uhh and we arrive at a relation between the 
upper left blocks of the two Lagrangians: 

Uhh = a^UNNia - bd'^c) = [a^ - c^d^y^b^UNNa. (53) 

This gives a relation between the operator forms of the Lagrangians that remain after in- 
tegrating out the respective lower components. To (9(l/m^) this expression is simplified 
to 

Uhh = a'UNNa. (54) 
The explicit expression as provided by Eq. (H5l) is 



Uhh =(l + ^BC] Unn U + ^Bc] , (55) 

which is the general form of the wave function renormalization to this order. It repro- 
duces, e.g., the special cases considered in [l^. Incidentally, this expression explains why 
the 0{l/m) expressions look equivalent — the transformation between the HF and FW La- 
grangians is (9(l/m^). 

The 1/m corrections to contact terms can be different for the two methods. This can be 
seen in, e.g., the ttN Lagrangian that has been developed to 0{l/m^) 19, 25] and higher [2^ 
in HBxPT, using the HF method. We find that the g\/Tn? corrections to the LECs di2 and 
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di3 in [19|, l20| are different in tlie HF and FW formalisms. Tlius, if these LECs are calculated, 
e.g., by resonance saturation, or extracted from data using the relativistic formulation, it 
would seem like the corresponding 1/m corrections would matter. This difference could be 
absorbed by a field renormalization and unitary transformation, see Refs jl^. I21I. provided 
the lower components of the relativistic Lagrangian are known. Therefore, the equivalence of 
the HF and FW approaches can be established quite generally and in particular regarding 
the values of LECs at higher orders. Since the FW result reproduces the relativistic S- 
matrix 17|, it is clear that so will the corresponding HF result. 



III. CONCLUSIONS 



We have carried out the heavy-fermion (HF) and Foldy-Wouthuysen (FW) non-relativistic 
reduction schemes to (9(l/m^), starting from a quite general, relativistic, Lagrangian. At 
this order explicit differences start to appear in the expressions of the respective non- 
relativistic Lagrangians. Our investigation treats a more general Lagrangian than the one 
in and is especially concerned with the LECs (introduced in the relativistic Lagrangian) 
as they appear in the NR theory, e.g., xPT. 

We have derived the field transformation between the two descriptions to (9(l/m^) and 
have shown how the 1/m corrections to certain LECs differ between the two approaches at 
this order. Although the Lagrangians differ at (9(l/m^), they can be related through a field 
redefinition, which we give in a more general form than the one in Ref. [l^. Thus, when 
using LECs of higher order, one must ascertain that the fields are appropriately redefined. 

We have also shown the equivalence between the NR propagators derived in HF and from 
an 1/m expansion. On the other hand, even though the HF and FW propagators go on-shell 
for the same four-momentum, they will differ off-shell. Also this difference is captured by 
the field redefinition. In order to do the field redefinition it is necessary to know the lower 
components of the original, relativistic Lagrangian. 
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APPENDIX: INTEGRATING OUT SMALL COMPONENTS 



We repeat here the arguments given in Ref. (see also Q) 

showing that the inte- 
gration over the lower components of the HF Lagrangian gives a constant. The discussion 
automatically carries over to the FW case. The path integral we need to perform for HF is 

JjjF = j[dh'^][dh']exp i j d^xh'\2m + iVo - ie)h' 

= Det(2m + ido + iTo - ie) = Det(2m + ido)I)et (l + h iTt 

= Det(2m + i^o) exp 



-Trln ( 1 + ^- iTo 

2 \ 2m + ioQ — ie 



2m + ido ~ 
= Const. 



(A.l) 



As argued in Ref. [3(], the propagator in the logarithm is always propagating in the backward 
time direction only. Thus the trace can not be closed, it has to vanish, and the determinant 
is a constant as indicated. 
To 0(l/m2) in FW we have 



Xpw = J [dn'^] [dn] exp 



i I d^xn^ {2m + iVo + A + -^{-^^ ■ ^ + C){-ia ■ V + B) 



1 



4m^ 
1 

8m 



{-la ■ V + C)(zPo + A){-'ia ■ V + B) 
^{{-ia ■ V + C){-ia -V + B), {iVo + D)}-ie]n 



(A.2) 



The resulting determinant can, just as in the HF case, be factorized into the free propagator 
and interaction terms: 

Xfw = Det I 2m + ido — — - + 



X Det 



1 + 



2m 8m^ 
1 



2^ + ^^o-£-^ 



[interaction terms) 



Det {2m + ido- 7^-7 + . . . 



2m Sm'^ 



X exp 



Tr In 1 + 



2'--'\' ■ 2m + ^^o-^-i^ + ...-^e 



'interaction terms) 



Const., 
(A.3) 



where in the last step we again use that the propagator propagates in the backward time 
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direction only, so that no closed loops are possible. 
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